Phase-space approach to the study of decoherence in quantum walks 
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We analyze the quantum walk on a cycle using discrete Wigner functions as a way to represent 
the states and the evolution of the walker. The method provides some insight on the nature of the 
interference effects that make quantum and classical walks different. We also study the behavior of 
the system when the quantum coin carried by the walker interacts with an environment. We show 
that for this system quantum coherence is robust for initially delocalized states of the walker. The 
use of phase-space representation enables us to develop an intuitive description of the nature of the 
decoherence process in this system. 



I. INTRODUCTION 

Quantum walks have been proposed as potentially 
useful components of quantum algorithms Q ■ In recent 
years these systems have been studied in detail and some 
progress has been made in developing new quantum algo- 
rithms using either continuous Q or discrete Q versions 
of quantum walks. The key to the potential success of 
quantum walks seems to rely on the ability of the quan- 
tum walker to efficiently spread over a graph (a network 
of sites) in a way that is much faster than any algorithm 
based on classical coin tosses. 

Quantum interference plays an important role in quan- 
tum walks being the crucial ingredient enabling a faster 
than classical spread. For this reason, some effort was 
made in recent years in trying to understand the impli- 
cations of the process of decoherence for quantum walks 
EM El 13 ■ Decoherence, an essential ingredient to under- 
stand the quantum-classical transition |8| , could turn the 
quantum walk into an algorithm as inefficient as its clas- 
sical counterpart. The models studied in this context can 
be divided in two classes depending on how the coupling 
with an external environment is introduced. In fact, a 
quantum walk consists of a quantum particle that can oc- 
cupy a discrete set of points on a lattice. In the discrete 
version, the walker carries a quantum coin, which in the 
simplest case can be taken as a spin- 1/2 degree of free- 
dom. The algorithm proceeds so that the walker moves 
in one of two possible directions depending on the state 
of the spin (for more complex regular arrays, a higher 
spin is required). So, in this context it is natural to 
consider some decoherence models where the spin is cou- 
pled to the environment and others where the position 
of the walker is directly coupled to external degrees of 
freedom. The specific system in which the algorithm is 
implemented in practice will dictate which of these two 
scenarios is more relevant. Several experimental propos- 
als to implement discrete quantum walks in systems such 
as ion traps , cavity QED ^} , and optical lattices [Tlj 
have been analyzed (see also Ref. 01 f° r a recent NMR 
implementation of a continuous quantum walk). 

The main effect of decoherence on quantum walks is 
rather intuitive: as the interaction with the environment 



washes out quantum interference effects, it restores some 
aspects of the classical behavior. For example, it has 
been shown that the spread of the decohered walker be- 
comes diffusion dominated proceeding slower than in the 
pure quantum case. This result was obtained both for 
models with decoherence in the coin and in the position 
of the walker [MIH0- However, it is known that classical 
correspondence in these systems has some surprising fea- 
tures. For example, for models with some decoherence in 
the quantum coin the asymptotic dispersion of the walker 
grows diffusively but with a rate that does not coincide 
with the classical one |(|. Also, a small amount of de- 
coherence seems to be useful to achieve a quantum walk 
with a significant speedup 0, Q ■ 

In this work we will revisit the quantum walk on a cycle 
(and on a line) considering models where the quantum 
coin interacts with an environment. The aim of our work 
is twofold. First we will use phase-space distributions 
(i.e., discrete Wigner functions) to represent the quan- 
tum state of the walker. The use of such distributions 
in the context of quantum computation has been pro- 
posed in Ref. 0] , where some general features about the 
behavior of quantum algorithms in phase space were no- 
ticed. A phase-space representation is natural in the case 
of quantum walks, where both position and momentum 
play a natural role. Our second goal is to study the true 
nature of the transition from quantum to classical in this 
kind of model. We will show that models where the en- 
vironment is coupled to the coin are not able to induce a 
complete transition to classicality. This is a consequence 
of the fact that the preferred observable selected by the 
environment is the momentum of the walker. This ob- 
servable, which is the generator of discrete translations in 
position, plays the role of the "pointer observable" of the 
system [8|,[l4|. Therefore, as we will see, the interaction 
with the environment being very efficient in suppressing 
interference between pointer states preserves the quan- 
tum interference between superpositions of eigenstates 
of the conjugate observable to momentum (i.e., position). 
Again, the use of phase-space representation of quantum 
states will be helpful in developing an intuitive picture 
of the effect of decoherence in this context. The paper 
is organized as follows: In Sec. II we review some ba- 
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sic aspects of the quantum walk on the cycle. We also 
introduce there the phase-space representation of quan- 
tum states for the quantum walk and discuss some of the 
main properties of the discrete Wigner functions for this 
system. In Sec. Ill we introduce a simple decoherence 
model and show the main consequences on the quantum 
walk algorithm. In Sec. IV we present a summary and 
our conclusions. 



II. QUANTUM WALKS AND THEIR 
PHASE-SPACE REPRESENTATION 

A. Quantum walks on the cycle 

The quantum walks on an infinite line or in a cycle with 
N sites are simple enough systems to be exactly solvable. 
For the infinite line the exact solution was presented in 
Ref . . The case of the cycle was first solved in Ref . . 
However, the exact expressions are involved enough to 
require numerical evaluation to study their main features. 
Here we will review the main properties of this system 
presenting them in a way which prepares the ground to 
use phase-space representation for quantum states (we 
will focus on the case of a cycle, the results for the line 
can be recovered from ours with t < N). 

For a quantum walk in a cycle of N sites, the Hilbert 
space is Tt = TLn®'H2, where Hat is the space of states of 
the walker (an ./V-dimensional Hilbert space) and H2 is 
the two-dimensional Hilbert space of the quantum coin 
(a spin 1/2). The algorithm is defined by a unitary evolu- 
tion operator which is the iteration of the following map: 

U walk = U°* H. (1) 

Here H is the Hadamard operator acting on the Hilbert 
space of the quantum coin [H = (o x + a z )j v2, Oi being 
the usual 2x2 Pauli matrices]. The operator U is the 
cyclic translation operator for the walker, which in the 
position basis is defined as U\n) = \n + 1) (mod N). It 
is worth stressing that the operator U a " is nothing but 
a spin-controlled translation acting as U a ~\n) <g> \ x ) = 
\n ± 1) <g> I ®). So, the map U wa ik consists of a spin- 
controlled translation preceded by a coin toss, which is 
implemented by the Hadamard operation (the use of the 
Hadamard operator in this context is not essential and 
can be replaced by almost any unitary operator on the 
coin Q). 

The notion of phase-space is natural in the context of 
this kind of quantum walk. In fact, the position eigen- 
states {\n), n = 0, . . . , N — 1} form a basis of the 
walkers' Hilbert space Hn- The conjugate basis is the 
so-called momentum basis {\k}, k = 0, TV — 1}. Po- 
sition and momentum bases are related by the discrete 
Fourier transform, i.e., 

(n\k) = -J=cxp{2ni?ik/N). (2) 



The cyclic translation operator U that plays a central 
role in the quantum walk is diagonal in momentum basis 
since U\k) — exp(— i2irk/N)\k). This simply indicates 
that momentum is nothing but the generator of finite 
translations. As a consequence of this, the total unitary 
operator defining the quantum walk algorithm is also di- 
agonal in such basis. This fact, which was noticed before 
by several authors, enables a simple exact solution of the 
quantum walk dynamics. Indeed, we can write the ini- 
tial state of the system using the momentum basis of the 
walker as (below p denotes the total density matrix of 
the system formed by the walker and the coin) 

iV-l 

p(0)= c k , k ,\k)(k'\®p 2 (0), (3) 

fe,fe'=0 

where P2(0) is the initial state of the quantum coin (we 
assume that the initial state of the combined system is a 
product, but this assumption can be relaxed). After t it- 
erations of the quantum map the reduced density matrix 
of the walker (denoted as p w ) is 

2V-1 

Pw (t)= c k ,vf(k,k',t) \k)(k'\, (4) 

fc,fc'=0 

/(fc,A;',t) = Tr 2 [M^ 2 (0)M fe t , t ] J (5) 

where the operator M k is defined as 

M k = exp(-i27r ka z /N) H. (6) 

All the temporal dependence is contained in the func- 
tion f(k,k',t) which can be exactly computed in a 
straightforward way: One should first expand the ini- 
tial state P2 (0) in a (nonorthogonal) basis of operators of 
the form |^1}(</>S| (l,m = 1,2), where \4> l k ) are the eigen- 
states of the operator M k (i.e., M k \(j) l k ) = \ k ,i\<j> l k ))- The 
explicit expressions for the eigenstates \<fi k ) and the eigen- 
values X k _i will not be given here since they can be found 
in the literature (see, for example, After doing this 

the evolution of the quantum state is fully determined by 
the equation 

f(k,k',t)= J2 (*w**,J' <4lMo)|<^> 

;,m=l,2 

(7) 

Below we will describe the properties of this solution us- 
ing a phase-space representation for the quantum state 
of the system. 

B. Phase-space representation 

Wigner functions are a powerful tool to represent 
the state and the evolution of a quantum system. For 
systems with a finite-dimensional Hilbert space, the dis- 
crete version of Wigner functions was introduced using 
different methods (see Ref. 01)- We will follow the 
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approach and notation used in Ref. |l3j, where these 
phase-space distributions were applied to study proper- 
ties of quantum algorithms. For completeness, we will 
give here the necessary definitions and outline some of 
the most remarkable properties of the discrete Wigner 
functions. 

For a system with an iV-dimcnsional Hilbert space the 
discrete Wigner function can be defined as the compo- 
nents of the density matrix in a basis of operators defined 
as 



A(q,p) = U q RV- p exp(inpq/N). 



(8) 



These are the so-called phase-space point operators. 
They arc defined in terms of the cyclic shift U (which 
in the position basis acts as U\n) = \n + 1)), the re- 
flection operator R (which in the position basis acts 
as R\n) = | — n)), and the momentum shift V (which 
generates cyclic displacements in the momentum basis, 
i.e., V\k) = \k + 1}). Phase-space operators are uni- 
tary, Hermitian and form a complete orthogonal basis 
of the space of operators (they are orthogonal in the 
Hilbert-Schmidt inner product since they satisfy that 
Tr[A(q,p)A(q' ,p')] = Nd qtq id PtP i). Expanding the quan- 
tum state in the A(q,p) basis as 



N-l 



P=Y1 W{q,p)A{q,p), 



(9) 



q,p=0 



the coefficients W(q,p) are the discrete Wigner functions 
of the quantum state, which are obtained as 
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W(q,p) = -Tr[pA(q,p)]. 



(10) 



This function has three remarkable properties that al- 
most give it the status of a probability distribution. The 
first two properties arc evident: Wigner functions are 
real numbers (a consequence of the Hcrmiticity of the 
phase-space operators) and they provide a complete de- 
scription of the quantum state (a consequence of the com- 
pleteness of the basis of such operators). The third prop- 
erty is less obvious: marginal probability distributions 
can be obtained by adding values of the Wigner func- 
tion along lines in phase-space. For this to be possible, 
it turns out that the phase-space has to be defined as a 
grid of 2N x 2N points where W(q,p) is given at each 
point precisely by Eq. JTHJ. Thus, adding the values of 
the Wigner function over all points satisfying the condi- 
tion aq — bp = c one gets the probability to detect an 
cigenstate of the operator D(b, a) = U V a exp(iirab/N) 
with eigenvalue exp(i7rc/2iV) (the sum is equal to zero 
if such eigenvalue does not exist). In particular, adding 
the Wigner function along vertical lines q — c one ob- 
tains the probability to detect eigenstates of the opera- 
tor D(0, 1) = V, with eigenvalues given by exp(i7rc/2iV). 
These numbers are equal to zero if c is odd and they are 
equal to the probability for measuring the position eigen- 
state \c/2) when c is even. Complementary, adding values 



of the Wigner function along horizontal lines enables us 
to compute the probability to detect a momentum eigen- 
state. 

A final remark about properties of the discrete Wigner 
function is in order. Figure ^ shows the Wigner function 
of a position eigenstate |no) and of a superposition of two 
position eigenstates, such as (|ni) + |ri2))/v / 2. As we see, 
in the first case the function is positive on a vertical line 
located at q = 2n and is oscillatory on a vertical line 
located at q = 2n + N. The interpretation of these os- 
cillations is clear. It is well known that Wigner functions 
display oscillatory regions whenever there is interference 
between two pieces of a wave packet. In this case, the 
cyclic boundary conditions we are imposing (that orig- 
inate from the fact that U and V are cyclic shift op- 
erators) generate a mirror image for every phase-space 
point. Thus, the oscillating strip can be interpreted as 
the interference between the positive strip and its mirror 
image. For the case of a quantum state which is a su- 
perposition of two position eigenstates, we observe two 
positive vertical lines with the usual interference fringes 
in between them. All these vertical lines have their cor- 
responding oscillatory counterparts originated from the 
boundary conditions which are located at a distance N. 
In what follows we will show Wigner functions for typical 
states of a quantum walker. 
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FIG. 1: Wigner functions W(q,p) for a localized state (up) 
and for a delocalized superposition state (down) . The dimen- 
sion of the Hilbert space is N = 41. The horizontal (vertical) 
axis corresponds to position (momentum). Color code is such 
that red (blue) regions correspond to positive (negative) val- 
ues while white corresponds to zero (no = 8, »n = 6, n,2 = 9). 
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For an initial state where the walker starts at 
a given position and the spin is initially unbiased 
[| V(0)> 2 = (|0) +i|l))/\/2], the behavior of the quantum 
walker starts to deviate from its classical counterpart at 
early times (in this paper we will only consider unbiased 
initial states for the quantum coin). The phase-space 
representation of the state is shown in Fig. [3| and makes 
evident that a peculiar pattern of quantum interference 
fringes develops between the different pieces of the wave 
packet. The consequence of these interference effect is 
evident also when one looks at the probability distribu- 
tion for different positions. This distribution is shown in 
Fig. Eland has been previously studied in the literature 
(see 0, 0, 0, E])- Like its classical counterpart, at a 
given time t the state initially located at hq has support 
only on states n satisfying that n + t + jiq adds up 
to an even number. However, in general the quantum 
distribution differs from the classical one, exhibiting 
peaks located at n = ±t/y/2 and a plateau of height 
l/V2t around no. After some time the Wigner function 
of the quantum walker develops a shape that resembles 
a thread, as it is clear in the pictures. For this reason 
we will call this a thread state. 
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FIG. 2: Probability distribution in position after t = 100 
iterations for an initially localized walker with unbiased spin. 
N = 301, no = 150. We only plot the function for sites such 
that n+t+no adds to an even number (solid), and also include 
the classical distribution (dotted). 



It is also interesting to analyze the evolution of the 
quantum walk for delocalized initial states. In particular, 
we will consider an initial state that is a coherent super- 
position of two position eigenstates (whose Wigner func- 
tion was already displayed in Fig. We find a Wigner 
function that develops into a sum of two threads with 
a region in between where interference fringes are evi- 
dent. This is displayed in Fig. 0] Some properties of the 
quantum walk for this kind of delocalized initial states 
were analyzed in Ref. where it was noticed that the 
asymptotic probability distribution can be rather differ- 
ent from the one obtained from a localized initial state. 
Below, we will show that the process of decoherence af- 
fects localized and delocalized initial states in a rather 
different way. 



III. DECOHERENCE AND THE TRANSITION 
FROM QUANTUM TO CLASSICAL 

A. The decoherence model 

We will consider a quantum walk where the quantum 
coin couples to an environment. To describe such cou- 
pling wc will use a model which was introduced and stud- 
ied in detail in Ref. [l^ . In that paper it was shown that 
one can mimic the coupling to an external environment 
by a sequence of random rotations applied to the quan- 
tum coin, which have the effect of scrambling the spin 
polarization. More precisely, these kicks will be gener- 
ated by the evolution operator 



exp(—iej n ■ a), 



(11) 



where the angles ej take random values and n is a fixed 
vector specifying the rotation axis. The virtue of this 
model is not only its simplicity but also the fact that can 
be experimentally implemented in a controllable manner 
using, for example, NMR techniques. 

After the application of one step of the quantum walk 
algorithm and one kick the evolution of the total system 
is 



p(t + l)= Kj U a ~~H p(t) HXJ- a ~~ K\, 



(12) 



To obtain a closed expression for the reduced density 
matrix of the walker for an ensemble of realizations of 
the random variables tj wc follow the method proposed in 
Ref. (see Ref. @ for a similar approach): Assuming 
that these angles are randomly distributed in an interval 
(— a,+a), this density matrix is 



r f de t r dc t 
p w {t) = Tr 2 / -^K t U^H...K 1 U^H 

J-a 2a J_ a 2a 

p w (0) <g> p 2 (0) HU- a -~K\ . . . HU-^Kl] (13) 

Expanding the initial state in the momentum basis as 
before enables us to simplify this expression. In fact, after 
doing this one can integrate over the random variables to 
find 



N-l 



Pw{t) 



J2 Ck,k'fn(k,k',t)\k}{k'\, (14) 
fc,fc'=0 



f n (k,k',t) = Tr 2 [O^(p 2 (0))]. 



(15) 



Here O n is a superoperator acting on the spin state (de- 
pending on direction n of the kicks) as 



O n ( P2 ) = ^-^-M kP2 Ml 



2 

(1-7) 



a n M k p 2 Ml,a. 



(16) 



where 7 = sin(2a)/(2a) is a parameter related to the 
strength of the noise (notice that 7=1 corresponds to 
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a) t = 2 b) t = 4 c) t = 8 

FIG. 3: Discrete Wigner function W(q,p) at different times. The initial state of the walker is localized at no = 20 (the dimension 
of the Hilbert space is N — 41). The quantum coin is in an unbiased initial state. Horizontal (vertical) axis corresponds to 
position (momentum). 
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a) t = 2 b) t = 3 c) t = 6 

FIG. 4: Discrete Wigner function W(q,p) at different times, for an initial state of the walker which is a superposition of two 
position eigenstates (N = 81, ni = 28, n% = 52) and an unbiased initial quantum coin. In (a) we plot the complete Wigner 
function (horizontal axis and vertical axis correspond to position and momentum, respectively). In (b) and (c) we only plot the 
Wigner function in the phase-space region defined by the black rectangle shown in (a). In this way the small scale oscillations 
of the Wigner function can be observed in detail. 



unitary evolution, i.e., to a = 0). One can find a sim- 
ple matrix representation for the superoperator O n for 
different choices of the rotation axis by writing p2(0) in 
the basis formed by the identity and the Pauli operators. 
In the Appendix we show the explicit form of this ma- 
trix representation, which is helpful in finding exact and 
numerical solutions to the problem. In what follows we 
will show results for the case h = y (the other cases are 
qualitatively similar). To find the state of the walker at 
arbitrary times we simply need to find eigenstates and 
eigenvalues of the superoperator O n . This can always be 
done numerically and also analytically in the interesting 
case of 7 = 0, which can be denoted as "total decoher- 
ence" . In such case, the exact solution turns out to be 

f y (k,k',t) = cos*(27r(fc -k')/N) 

/ sin( 2 ^-fc')/AQ \ (f) 



Several features of the decoherence effect are evident in 
the above formula. The environment produces a ten- 
dency towards diagonalization of the density matrix of 
the walker in the momentum basis (matrix elements with 
k — k' = N 1 4 are maximally suppressed) . The decay of 
nondiagonal elements is exponential in time, as already 
discussed in Ref. It i s a l so clear that momentum 

eigenstates are not affected by the interaction since they 
are eigenstates of the full evolution [in fact, from Eqs. 
(I15fl and (|l(j|) follows f n (k,k,t) = 1]. In this sense they 
are perfect pointer states for this model. In what follows 
we will present results concerning the evolution of several 
initial quantum states. 
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FIG. 5: Discrete Wigner function W(q,p) for a fixed time (t = 11) and different values of the strength of the coupling to the 
environment (a). The initial spin is in an unbiased state. Horizontal and vertical axis respectively correspond to position and 
momentum (total dimension of the Hilbert space is N = 41 and the initial state is located at no = 20). 




FIG. 6: Probability distribution in position at t = 100, for different values of the coupling to the environment (parametrized 
by a). The initial state of the quantum coin is unbiased (iV = 301, no = 150, and n = y). We only plot the function for sites 
such that n + t + no adds to an even number (solid), and also include the classical distribution (dotted). Fig. |5|shows the same 
plot without decoherence. 



B. Decoherence, pointer states, and the transition 
from quantum to classical 

The effect of decoherence on the evolution of states 
which arc initially localized in position has been ana- 
lyzed elsewhere 0,1a, • As shown in Fig. the Wigner 
function of the evolved quantum state gradually loses its 
oscillatory nature. Thus, instead of a thread state the 
interaction with the environment gradually produces a 
mixed state with a binomial distribution in the position 
direction (which has an approximately Gaussian shape 
for large t) but remains constant along the momentum 
direction. It is worth noticing that for any value of 7 the 
resulting state has support only on position eigenstates 
satisfying that the sum n + n + t is equal to an even 
number, as it was already pointed out for both the clas- 
sical distribution (7 = 0) and the purely quantum one 
(7 = 1) in the preceding section. It is interesting to no- 
tice that the process of decoherence has a rather simple 
interpretation when represented in phase-space: Deco- 



herence in phase-space is roughly equivalent to diffusion 
in the position direction. This is not unexpected: In fact, 
in ordinary quantum Brownian motion models a coupling 
to the environment through position (momentum) gives 
rise to a momentum (position) diffusion term in the evo- 
lution equation for the Wigner function. The situation 
here is quite similar, since the walker effectively couples 
with the environment through its momentum. This is 
indeed the case because the environment interacts with 
the quantum coin which is itself coupled with the walker 
through the displacement operator which is diagonal in 
the momentum basis. Therefore, the decoherence effect 
on the Wigner function is expected to correspond to dif- 
fusion along the position direction. 

As noticed before, if one considers initial states where 
the walker is localized in a well defined position, one 
can see that the probability distribution for the different 
positions of the walker gradually tends to the classical 
one by increasing the coupling strength a from a — 
(no kicks) to a = tt/2. This is shown in Fig. 
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FIG. 7: Absolute value of the density matrix in the position basis (left) and discrete Wigner function W(q,p) (center) for a 
state evolving from an initially delocalized state of the quantum walker under full decoherence (t = 6, a = 0.5tt, ni = 38, 
ri2 = 62, N = 101). The presence of quantum interference is manifested in the nondiagonal terms of the density matrix and in 
the oscillations of the Wigner function. The small scale oscillations of the Wigner function are shown in the right plot, which 
shows W(q,p) in the phase-space region defined by the black rectangle depicted in the center plot. 



From the above analysis we could be tempted to con- 
clude that the interaction between the quantum coin 
and the environment induces the classicalization of the 
walker. However, this is not the case. The process of 
decoherence induced in this way is not complete. This is 
most clearly seen by analyzing how it is that the inter- 
action with the environment affects initial states of the 
quantum walker which are not initially localized. In Fig. 
0we show the Wigner function of an initially delocalized 
state (shown in Fig. ^) under full decoherence. We can 
clearly see that decoherence does not erase all quantum 
interference effects. In fact, as mentioned above, the in- 
teraction with the environment induces diffusion along 
the position direction. Therefore, interference fringes 
which arc aligned along the position direction are im- 
mune to decoherence. Thus, the final state one obtains 
from a superposition of two position eigenstates is not the 
mixture of two binomial states but a coherent superposi- 
tion of them. This peculiar behavior is easily understood 
by noticing that this is a simple consequence of the fact 
that momentum eigenstates are pointer states: Decoher- 
ence is effective in destroying superpositions of pointer 
states but highly inefficient in destroying superpositions 
of eigenstates of the conjugate observable (position). 

C. Entropy 

By analyzing the entropy of the reduced density ma- 
trix of the walker one can get a more quantitative mea- 
sure of the degree of decoherence achieved as a conse- 
quence of the interaction with the environment. For 
convenience we will not examine the von Neumann en- 
tropy Sy but concentrate on the linear entropy defined 
as Sl = ~ m (Tr[p^]), which is easier to calculate. This 
entropy provides a lower bound to Sy |20j . It is possible 
to show that almost no entropy is produced by the decay 



of the coherence present in the initially delocalized su- 
perposition state. In fact, this can be seen by comparing 
the entropy produced from the initially delocalized su- 
perposition and the one originated from an initial state 
in which the walker is prepared in an equally weighted 
mixture of two positions. These entropies can be seen 
in Fig. [5] The initial entropy of the mixture is 1 bit 
[ln(2)]. It is quite clear from the curves shown in such 
figure that the entropy arising from the initial mixture 
remains to be 1 bit higher than the one originated from 
the initial coherent superposition. Thus, the quantum 
coherence present in the initial state is robust under the 
interaction with the environment and does not decay at 
all. 

Fig. |S| shows another interesting feature: One would 
naively expect a monotonic dependence of the en- 
tropy with the coupling to the environment (which is 
parametrized by a). However this is not the case since 
the curves in Fig. |S| intersect. This peculiar effect is 
made more evident in Fig. [5] where we study the entropy 
at a fixed time as a function of the coupling strength. 
In this figure a clear indication of a nontrivial behavior 
is seen: For early times the entropy grows slowly with a 
and exhibits a flat plateau for large values of a. However, 
as time progresses a peak develops: The largest value of 
entropy at a given time is not achieved by the largest 
coupling. To the contrary, the largest entropy is attained 
by an intermediate coupling a c , whose value decreases 
with time. 

The fact that for a given time the maximal entropy 
is not achieved by the maximal coupling to the environ- 
ment is counterintuitive. As entropy is a measure of the 
spread of a distribution, this strange behavior can be 
rephrased as a manifestation of the counterintuitive fact 
that the decohered state (which is approximately diago- 
nal in position basis) has a probability distribution that 
is more spread for a = a c than for a > a c . A possible 
explanation for this peculiar behavior is the following: 
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For high values of the coupling to the environment the 
state rapidly becomes classical and the spread in position 
grows diffusively, as in the classical random walk. When 
the coupling to the environment is not strong, our result 
seems to indicate that the state of the walker remains 
"quantum" for a longer time during which it spreads at 
a rate faster than classical. When this quantum state 
finally decoheres it may end up having a larger entropy 
than the one attained for high coupling simply because 
it is spread over a wider range of positions. We specu- 
late that there could be a relation between this peculiar 
feature and the properties that make some degree of deco- 
herence useful for quantum walks as discussed by Kendon 
and Tregenna in 0, • The value of the a c introduced 
above depends on both N and t and could be related to 
the position of the minima reported in Refs. 0,0- For 
example, in a cycle regime (t » N) a c diminishes with 
increasing N as it is also the case for the position of the 
minima of the so-called quantum mixing time 0, 0, Q • 
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FIG. 8: Linear entropy as a function of time for various values 
of the system-environment coupling strength: a = (solid), 
O.Itt (dotted), 0.2tt (dash dotted), 0.5vr (dashed). Fhe top 
(down) plot corresponds to an initial state which is an equally 
weighted mixture (superposition) of two position eigenstates. 
N = 401, ni = 150, n 2 = 250. 



IV. CONCLUSIONS 

The use of phase-space representation enables us to 
develop some intuition about the nature of the decohcr- 
ence process in the kind of quantum walk analyzed in 




FIG. 9: Linear entropy as a function of the coupling strength 
a for various values of time: t = 5 (solid), 10 (dashed), 50 
(dotted), 100 (□), 300 (x), 500 (dash dotted). The initial 
state of the walker is well localized and the initial state of 
the quantum coin is unbiased. N = 401 and all the curves 
are below the saturation regime [ln(JV) w 5.994]. It is evident 
that, after some time, the maximum value of the entropy is 
not achieved by the maximum value of the coupling strength. 



this paper. By coupling the quantum coin to an environ- 
ment we obtain a decoherence model which is roughly 
equivalent to position diffusion. As we mentioned above, 
this is a natural result whose origin can be traced back 
to the way in which the system effectively couples to the 
environment (via the momentum operator). The relation 
between decoherence and position diffusion can also be 
established by analyzing in more detail the structure of 
the superoperator O n [given in Eq. I|lb|) ]. Let us consider 
the form of the superoperator after t iterations. If we use 
Eq. i|lb|) we can easily see that, as each iteration doubles 
the number of terms, we will have an expression with 
2* terms each one of which has Pauli operators applied 
at different times. To obtain the function f n (k, k',t) one 
should compute the trace over the quantum coin. In each 
of the 2* terms we can move the Pauli operator a n to- 
wards the outside of the expression and cancel them due 
to the cyclic property of the trace. For the case h = y 
it is easy to show that the only remaining effect of the 
Pauli operators (that in this case anti-commute with the 
Hadamard operator) is to reverse the direction of the ro- 
tation in Mk defined in Eq. © . The final expression can 
be shown to be 



Tr(0») 



E ••• E ^ a +?)*-* a -7)* 



Q t =0,l 

x Tr 



(M kt ...M klP Mi,...Ml^ (18) 

where ctj = a± + ... + ctj and kj = (— l) Qj ' _1 fc (we use 
the convention ao =0). Therefore, the superoperator is 
the sum of 2* terms each one of which contains a contri- 
bution that is identical to that of a quantum walk where 
the direction of the walker is chosen at random after the 
first step. Each of the 2* terms is labeled by a bit string 
(ai, . . . , a n ) and corresponds to a quantum walk where 
the direction of the jth step (j > 2) is reversed if and only 
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if ctj-i is odd. In the limit of total decoherence each of 
these terms has equal weight. Therefore the final state is 
simply the average over an ensemble where each member 
corresponds to each of the 2' -1 possible choices of two di- 
rections (forward or backward) for the t—1 steps (notice 
that the direction of the first step is not affected by the 
decoherence model we chose). For this type of decoher- 
ence it is clear that the quantum walk becomes a random 
walk. The relation between decoherence and position dif- 
fusion is quite evident in this way. It is worth pointing 
out that similar models of decoherence were considered 
in Ref. [2lJ in a different context. 

Other decoherence models have been analyzed for 
quantum walks 0, Q, El, where the effective coupling 
to the environment is through the position observable. 
In such case, we expect decoherence to correspond to 
diffusion along the momentum direction. Combining the 
two types of decoherence (i.e., considering coupling to 
the environment via the quantum coin and the position 
of the quantum walker) the initial state corresponding to 
a superposition of two positions would finally decay into 
a mixture of two binomial states (see Ref. [21| for similar 
results obtained when studying decoherence models with 
a natural phase-space representation in a finite quantum 
system evolving under various quantum maps). 

The above conclusions are generic for any model in 
which decoherence is due to the coupling of the quantum 
coin to an environment. An interesting class of models, 
based on the use of quantum multi-Baker maps |22j , has 
been studied. In such models one replaces the quantum 
coin with a quantum system with a higher-dimensional 
Hilbcrt space. The total space of states is then Tt = H.n® 
TLm- Here M is the dimensionality of the system which 
plays the role of the quantum coin and is considered to be 
an even number (M = 2m, so that we can always consider 
Hm = <£> H m ). The dynamics for a quantum multi- 
Baker map is defined in terms of the unitary operator 
[that replaces Eq. JIJ]: 

Umulti-Baker = U " Bm, (19) 

where Bm is the unitary operator defin ing the so-called 
"quantum Baker map" (see Refs. 0, |23|) and a z is a 
Pauli operator acting on the Hilbert space H.2 (the most 
significant qubit of the internal space Hm = 7^2 <8> 7~C m )- 
The properties of the o per ator Bm have been widely 
studied in the literature 23]: The map faithfully repre- 
sents a classically chaotic system (in the large M limit). 
From the point of view of the quantum walker the sit- 
uation is quite similar to the one we studied in this 
paper. One can describe a quantum multi-Baker map 
as an ordinary quantum walk where the quantum coin 
(whose Hilbert space is H.2) interacts with an environ- 
ment (whose Hilbert space is Ti m ). The interaction is 
modeled by the quantum Baker map acting on the to- 
tal internal Hilbert space %, which also replaces the 
usual Hadamard step in Eq. IjTjl. As the quantum Baker 
map is chaotic, the state of the quantum coin will be 
roughly randomized after each iteration. Thus, the effect 
should be similar to the one we described here (where 
the quantum coin is subject to a noisy evolution). How- 



ever, after a large number of iterations (of the order of 
M) all the possible orthogonal directions available in the 
internal space of the quantum coin would have been ex- 
plored. One should therefore expect that this model will 
stop being effective in producing decoherence after such 
time. Recent studies of quantum multi-Baker systems 
agree with these expectations (see Ref. where a tran- 
sition from diffusive to ballistic behavior after a time of 
the order of M has been analyzed). In any case, based 
on the results of our work we believe that in quantum 
multi-Baker systems the relative stability of initially de- 
localized superpositions will also be observable. 

We thank Marcos Saraceno and Augusto Roncaglia 
for useful discussions and assistance during several 
stages of this work. C.C.L. was partially supported 
by funds from Fundacion Antorchas, Anpcyt and Ubacyt. 



APPENDIX A: DENSITY MATRIX UNDER 
DECOHERENCE 

A matrix representation for the superoperator O n can 
be obtained for n = x,y, z. We will write this superop- 
erator in the basis formed by the identity and the Pauli 
matrices (we use the standard ordering of the basis as 
{/, a x , a y , <J Z }). In such basis the matrix of O n is 



On 



cos(— ) — isin(— ) 

7sin(+) 7COs(+) 

-cos(+) sin(+) 

— 7isin(— ) 7COs(— ) 



where [±) = (27r(fc ± k')/N). 

Similar results arc obtained for fi = x, z: 



O x 




o z = 



-isin(-) 

sin(+) cos(+) 

— 7Cos(+) 7sin(+ 

7cos(-) 



( cos(-) -isin(-) 

7sin(+) 7COs(+) 

— 7Cos(+) 7sin(+) 

\ — isin(— ) cos(— ) 



Notice they all converge to the same matrix when 7 = 1 
(no decoherence). 

To compute Eq. I|14|l we need to diagonalizc O n . Al- 
though the matrix representation of the superoperator is 
rather sparse, the eigenvalues and eigenvectors are quite 
cumbersome for an arbitrary value of 7 (including 7 = 1, 
no decoherence) , so it was more convenient to use numer- 
ical techniques. However, for the special case of complete 
decoherence (7 = 0) it is possible to obtain a simple for- 
mula and the final result for f n (k,k',i). The result for 
h = y, z is given in Eq. (|17|) . 
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